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ABSTRACT 

We study the Spectral Energy Distribution (SED) and the power spectrum of Galactic cirrus 
emission observed in the 14 deg 2 Science Demonstration Phase field of the He rschel-ACTLAS 
using Herschel and IRAS data from 100 to 500 fim. We compare the SPIRE 250, 350 and 
500 ^m maps with IRAS 100 ^m emission, binned in 6' pixels. We assume a modified black- 
body SED with dust emissivity parameter /3 (F oc A~P) and a single dust temperature Tj, and 
find that the dust temperature and emissivity index varies over the science demonstration field 
as 10 < Ta < 25 K and 1 < /3 < 4. The latter values are somewhat higher than the range of f3 
often quoted in the literature (1 < (3 < 2). We estimate the mean values of these parameters 
to be 7d = 19.0 ± 2.4 K and f3 — 1.4 + 0.4. In regions of bright cirrus emission, we find that 
the dust has similar temperatures with 7a = 18.0 + 2.5 K, and similar values of ft, ranging 
from 1.4 ± 0.5 to 1 .9 + 0.5. We show that and [3 associated with diffuse cirrus emission are 
anti-correlated and can be described by the relationship: /3(T^) = NT^ with [N — 116 ± 38, 
a = -1.4 + 0.1]. The strong correlation found in this analysis is not just limited to high 
density clumps of cirrus emission as seen in previous studies, but is also seen in diffuse cirrus 
in low density regions. To provide an independent measure of T& and /3, we obtain the angular 
power spectrum of the cirrus emission in the IRAS and SPIRE maps, which is consistent with 
a power spectrum of the form P(k) = Po(k/ko) y where y = 2.6 + 0.2 for scales of 50 - 200' in 
© 0000 RAS tne spire maps. The cirrus rms fluctuation amplitude at angular scales of 100' is consistent 

with a modified blackbody SED with 7^ = 20.1 ± 0.9 K and (3 — 1.3 ± 0.2, in agreement with 
the values obtained above. 
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1 INTRODUCTION 

The sub-millimeter and millimeter emission of diffuse Galactic 
dust is primarily determined by the thermal radiation of large dust 
grains that are in equilibrium with the interstellar radiation field 
(Desert et al 1990). The dust organizes itself into large-scale struc- 
tures such as cirrus and filaments both at low and high Galactic 
latitudes (Low 1984). To compare with previous studies we approx- 
imate the dust emission by a single thermal Planck spectrum at the 
temperature r d of the grains modified by a power law dependence 
on frequency parameterized by the spectral emissivity parameter p 
in the optically thin approximation: 

7(v) = e(v)B(v,r d )iV H , (1) 

where I(v) is the specific brightness, B Y is the Planck spectrum, 
and Njj is the total hydrogen column density along the line of sight. 
The isothermal assumption is likely to be only approximate. While 
at high latitudes the overlap of several dust sources along the line of 
sight is expected to be small, we still expect a complex temperature 
structure due to variations in the grain size distribution and any 
variations in the radiation field. In the above equation, e(v) is the 
emissivity 

6(y) = X rf e (^) (2) 

where X d is the dust-to-gas mass ratio and e is the emissivity at 
frequency v . 

With Ti at the level of 10 to 30K the maximum intensity is 
found at far-infrared/sub-mm wavelengths. Due to the lack of cov- 
erage at far-IR wavelengths, studies on the Galactic cirrus temper- 
ature and emissivity before Herschel (Pilbratt et al. 2010) focused 
on the wavelength bands shorter than 160 that were covered 
by IRAS and Spitzer (e.g., Miville-Deschenes et al. 2002; Jeong et 
al. 2005), longer than 1 mm covered by various cosmic microwave 
background (CMB) experiments (e.g., Desert et al. 2009; Veneziani 
et al. 2010), and a combination including limited sub-mm data (e.g., 
Dupac et al. 2003; Paradis et al. 2009; Bernard et al. 2009). 

With its unprecedented angular resolution and the ability to 
cover wide fields by scanning across the sky, the Spectral and Pho- 
tometric Imaging Receiver (SPIRE; Griffin et al. 2010) on Herschel 
now allows for the first time the possibility to study the structure of 
the ISM at tens of arcseconds to degree angular scales at the peak of 
the dust spectral energy distribution (SED). At the smallest angu- 
lar scales probed by SPIRE the structure of the diffuse interstellar 
medium provides information about the initial conditions for the 
formation of dense molecular clouds. Some of these clouds may be 
on the verge of gravitational collapse leading to the formation of a 
new star (Olmi et al. 2010, Sadavoy et al. 2010, Ward-Thompson 
et al. 2010). At large angular scales, dust is a key tracer of the 
large-scale physical processes occurring in the diffuse interstellar 
medium (Miville-Deschenes et al. 2010). Dust emission is also re- 
lated to the density structure of diffuse clouds and could potentially 
provide a way to study the projected density distribution within the 
Galactic cirrus. 

In equation (2), the dust emissivity index /3 provides informa- 
tion on the physical nature of dust and connects the grain struc- 
ture to the large-scale environmental density. The spectral index of 
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by European-led Principal Investigator consortia and with important partic- 
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the emissivity depends on grain composition, temperature distribu- 
tion of tunneling states, and the wavelength-dependent excitation 
(e.g., Meny et al. 2007). The emissivity f3 is also expected to vary 
with wavelength when the dust temperature is intrinsically multi- 
component but described by an isothermal model (Paradis et al. 
2009). While the dust SED models in the literature generally as- 
sume a fixed value for the spectral emissivity /? between 1.5 and 
2.5, there might be significant variations in ji, even in a small cirrus 
region, when taking into account the disordered structure of dust 
grains. There are also however the uncertainties associated with the 
dust size distribution and the silicate versus graphite fractions; large 
variations in both these quantities could result in different equilib- 
rium temperatures even in a small cirrus region. This complicates 
any physical interpretation of /? and 7d when using an isothermal 
SED. Instead of multi-component models, we use the isothermal 
model here so that we can compare our results to previous anal- 
yses that make the same assumption. Our SED modeling is also 
limited to four data points. Of particular interest to this work is 
the suggestion that /? and T& are inversely related in high density 
environments of Galactic dust by previous observations in the sub- 
millimeter and millimeter domain, both at low (Dupac et al 2003, 
Desert et al 2009) and high (Veneziani et al 2010) Galactic lati- 
tudes. 



Most of the recent studies on properties of the Galactic cir- 
rus focused on high density environments, such as cold clumps and 
molecular clouds with intensities of order 100 MIy sr~' or more at 
far-IR wavelengths (see, however, Bot et al. 2009 for a study on dif- 
fuse medium at small scales). Properties of the interstellar medium 
at high latitudes, especially involving diffuse cirrus with intensities 
of order a few MIy sr~', are still not well known. A good modeling 
of diffuse dust distribution and its characteristics in the high lati- 
tude regions is necessary in order to remove its contamination from 
CMB anisotropy measurements (see for example Leach et al. 2008 
and Ricciardi et al. 2010). The CMB community primarily relies 
on models developed with IRAS and DIRBE maps to describe the 
dust distribution (e.g., Schlegel et al. 1998) and the frequency de- 
pendence of the intensity (e.g., Finkbeiner et al. 1999). With wide- 
field Herschel-ATLAS (H-ATLAS; Eales et al. 2010) maps we can 
study the dust temperature and emissivity variation across large ar- 
eas on the sky, first in the SDP field covering 14 deg 2 and eventually 
over 550 deg 2 spread over 5 fields with varying Galactic latitudes. 
The H-ATLAS SDP patch is at a Galactic latitude of ~ 30 degrees. 
Combining SPIRE data at 250, 350, and 500 yum with IRAS maps 
of the same area at 100 /im allows us to sample the peak of the dust 
SED accurately. 



Here we present an analysis of diffuse Galactic cirrus in the 
H-ATLAS SDP field from 100 yum to 500 fim using IRAS and 
Herschel-SVIRE maps. We derive physical parameters of the dif- 
fuse dust at arcminute angular scales such as the temperature and 
spectral emissivity parameter and their relationship to each other. 
We also present a power spectrum analysis of the cirrus emission, 
which allows us to study the spatial structure of the interstellar 
medium from tens of arcseconds to degree angular scales. The dis- 
cussion is organized as follows: Section 2 describes the datasets; 
Section 3 describes the pipeline adopted. Section 4 reports results 
related to the temperature and spectral emissivity parameter while 
Section 5 describes the cirrus power spectrum. We conclude with a 
summary in Section 6. 
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Figure 1. Selected regions of high cirrus intensity that we have individually analyzed (labeled 1 to 5) as well as the two low-cirrus regions (labeled Fbl and 
Fb2) that we used to account for the extragalactic background intensity. The map, at 250 (im, is color coded in intensity units of Jy/sr. 



2 DATA SETS 

We use Hersc-hel-SPlRE maps in the H-ATLAS 14 deg. 2 Sci- 
ence Demonstration Phase (SDP) field, centered at RA=9 h 25 m 31 m , 
DEC=0 d 29 m 58 s , overlapping with the GAMA survey (Driver et al. 
2009). In addition to the three SPIRE bands, we also use the IRAS 
100 /im map. The latter is obtained by IDL routines projecting the 
HealpijQ format map made by the IRIS processing system of the 
IRAS surve)0 (Miville-Deschenes & Lagache 2005). 

We refer the reader to Pascale et al. (2010) for details on the H- 
ATLAS SPIRE map making procedure and basic details related to 
the maps. Since we are interested in the diffuse emission, we make 
use of a set of maps that have been especially made to preserve 
the extended structure by accounting for the map-making transfer 
function. We also produced a second set of maps using the same 
timelines processed by HIPE (Ott et al. 2006), but with an inde- 
pendent map-making pipeline. This involved the use of an itera- 
tive approach to make new maps using SHIM vl.O (The SPIRE- 
HerMES Iterative Mapper; Levenson et al. 2010). For that map 
maker simulations show a transfer function that is close to unity 
over arcminute to degree angular scales. The results we describe 



http://healpix.jpl.nasa.gov 
http://www.cita.utoronto.ca/ mamd/IRIS/ 



here, however, are consistent within overall uncertainties between 
the two sets of maps. Thus, we describe results primarily using the 
H-ATLAS map making pipeline of Pascale et al. (2010). We do 
not use the H-ATLAS PACS SDP maps for this analysis since the 
diffuse emission at short wavelengths imaged by PACS is heavily 
filtered out during the map-making process as carried out in the pro- 
duction of H-ATLAS PACS SDP data products (Ibar et al. 2010). 



3 DATA ANALYSIS 

The aim of this paper is to characterize the physical properties of 
dust emission over the whole SDP area. In order to avoid contami- 
nation of our Galactic dust measurements from extragalactic point 
sources, we first remove the bright detected sources from each of 
the maps making use of the H-ATLAS source catalogs (Rigby et al. 
2010). This catalog involves sources that have been detected at 5cr 
in at least one of the bands. Given that SPIRE data have beam sizes 
of approximately 18, 25, and 36", respectively, at 250, 350 and 500 
fim, we introduce a source mask by simply setting the pixel values 
over a square size of 20, 30, and 40" to be zero at the source loca- 
tions. For a handful of extended sources in the Rigby et al. (2010) 
catalog, we increased the size of the mask based on the source size 
as estimated directly from maps. With close to 6700 sources in to- 
tal, this masking procedure involved a removal of 2, 5 and 6% of the 
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Table 1. Dust intensity and SED properties of selected high-intensity features from the SDP map 



Region # 




i 


2 


3 


4 


5 


RA (deg) 




137.44 


136.18 


135.69 


134.79 


136.35 


DEC (deg) 




2.39 


1.90 


1.63 


1.47 


-0.31 


Area fripp ") 




0.33 


0.063 


0.29 


0.18 


0.14 


Flux 100 iim (MJy sr" 


') 


1.2±0.1 


1.5+0.1 


1.6±0.1 


1.2+0.1 


0.9±0.1 


Hux250//m(MJy sr" 


') 


2.4±0.4 


3.5±0.5 


3.3±0.5 


2.7±0.4 


1.8+0.3 


Mux 3MJ /im (MJy sr 


) 


1.1±U.Z 


1.8±U.3 


1.6±U.3 


1 .3+0. z 


n ixfl 1 
U. /±U. 1 


Flux 500 fim (MJy sr" 


') 


0.5±0.1 


0.9±0.1 


0.8±0.1 


0.7±0.1 


0.4±0.1 


In A 




-9.1±0.9 


-9.2±0.8 


-9.0±0.9 


-9.5±0.9 


-9.3±1.1 






1.8+0.5 


1.4±0.4 


1.6±0.5 


1.4±0.5 


1.9 ±0.6 


r d 




17.6±2.3 


18.3±2.2 


18.1+2.5 


18.3±2.5 


17.4±2.9 


* 2 




0.8 


0.5 


0.6 


1.1 


2.7 



data at 250, 350 and 500 yum, respectively. Removal of such a small 
fraction of pixels does not bias our results. Alternatively, we could 
have modeled each source by fitting the PSF at each of the source 
locations and removing the flux associated with the source and re- 
taining the background; we ran a set of simulations to study if the 
two approaches lead to different results, but we did not find any. In 
the case of the 100 fim IRAS map, we similarly masked roughly 35 
point sources in the SDP area from the IRAS-Faint Source Catalog 
of Wang & Rowan-Robinson (2009) by setting the pixel intensity 
to be zero over a square area of size 6'. We also found consistent 
results when we replace all IRAS pixels above 5 tr with zero inten- 
sity. 

To compare the IRAS and Herschel maps, we then convolve 
the source-masked SPIRE maps to the angular resolution of the 
IRAS 100 fim map (258" FWHM) and repixelize SPIRE maps at 
the same pixel scale as IRAS (120" pixel sizes). For the IRAS map, 
the flux errors are estimated by assuming that noise is isotropic 
and using the IRIS noise estimate (Miville-Deschenes & Lagache 
2005). For SPIRE, we use the noise maps produced by taking the 
differences of repeated scans (Pascale et al. 2010) and convolve the 
noise map to the IRAS 100 fim angular resolution and IRAS pixel 
size. 



4 DUST SPECTRAL ENERGY DISTRIBUTION 

To describe the SED we use a modified black-body spectrum with 
a spectral emissivity fl in the optically-thin limit such that: 

S(v)=A^ B v (T d ) (3) 

where the amplitude A depends on the optical depth through the 
dust, fi is the spectral emissivity, and r d is the temperature of the 
dust. We take vo = 3 THz corresponding to the IRAS 100 /im mea- 
surement. 

To estimate the best-fit values for the three unknown param- 
eters (A,/?, r d ), we make use of a Markov Chain Monte Carlo 
(MCMC) approach (Lewis & Bridle 2002). We are able to gen- 
erate the MCMC chains rapidly and at the same time fully sample 
the likelihood functions of the model parameters. Appropriate sam- 
pling of the likelihood is crucial to study the relation between /3 and 
T d (Section 4.4). 

The model fits adopt uniform priors over a wide range with 
P between -1 and 6 and r d between 1 and 50 K so that we do not 



incorrectly constrain the best-fit values by a narrow range of priors. 
The ranges are set such that we also allow both the minimum and 
maximum values of fi and T d to be well outside the expected ex- 
tremes. We run MCMC chains for each region until we get conver- 
gence based on the Gelman and Rubin statistic (Gelman & Rubin 
1992), with a value for 1 - R of at least 0.01 where R is defined as 
the ratio between the variance of chain means and the mean of the 
variances. 

4.1 Zero-level in SPIRE maps 

Since the SPIRE maps are not absolutely calibrated, to study the 
dust temperature over the SDP region as a whole, we account 
for the zero-point offset through the pixel-correlation method of 
Miville-Deschenes et al. (2010). Before computing the correlation 
we remove a constant intensity corresponding to the extragalactic 
background at 100 /jm from the IRAS map (0.78 MJy sr~'; Lagache 
et al. 2000; Miville-Deschenes et al. 2007). We then correlate the 
IRAS pixel intensity with that of a SPIRE map at the corresponding 
pixel. Here we use the source-masked SPIRE maps repixelized to 
the IRAS pixel scale following the procedure described in Section 3. 

We show this correlation in Figure 2 where we plot the pixel 
intensity values at the three SPIRE bands as a function of the IRAS 
intensity. We minimize the difference between Hy[Sy(/l) - Py(/l)] 
where Sy(/1) is the SPIRE intensity at pixel (i,f) at waveband A 
and Pjj(A) is the predicted SPIRE flux in each of the pixels by scal- 
ing IRAS 100 nm map intensity /(100), with the assumption that 
PijW = G x /y(100) + So, where G, sub-mm color (also called 
"gain" in Miville-Deschenes et al. 2010), and So, the zero-point 
offset, are the mean values over the whole of the map. This correct 
for the additive (So) term under the assumption that the IRAS map is 
a true reflection of the sky. The estimated values of G are 1.9 + 0.3, 
0.9 ± 0.3, and 0.5 ± 0.1 at 250, 350, and 500 /urn respectively, while 
S takes the values of -(3.2 ± 1.0), -(1.5 ± 0.6), -(0.8 ± 0.4) MJy 
sr~' at each of the three frequencies. 

The uncertainty in sub-mm color G is estimated by taking the 
rms of the ratio involving S (A)/I( 100) once the offset S is removed 
from each of the maps, while the uncertainty in So comes from the 
rms of the difference involving S(A) - Gx /(100), with the best-fit 
values for sub-mm color used in the computation. As described in 
Miville-Deschenes et al. (2010) these rms values reflect the overall 
uncertainties more accurately than if one were to simply use the 
statistics associated with the linear fit to the relation between S (A) 
and /(100). Those have uncertainties that are at least a factor of 10 
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Figure 2. The /RAS 100 /mi intensity vs. SPIRE intensity at 250 (top), 350 
(middle), and 500 (bottom) fim. The SPIRE maps are pixelized to IRAS and 
we use the relation between SPIRE intensity and IRAS intensity to deter- 
mine the sub-mm color and zero-point offset relative to the IRAS map (see 
text for details). The points colored blue are from the bright cirrus regions 
marked 1 to 5 in Figure 1. The points with zero values to SPIRE inten- 
sity are pixels that remain zero after smoothing to IRAS pixel scale due to 
our point and extended source masks that were defined in the original pixel 
scale. These and other points that fall off the main band of points, possi- 
bly associated with issues related to parts of the time-streams that are either 
contaminated or contain glitches that were removed, do not bias the color 
and zero-point offset estimates. 



Figure 3. The Galactic dust SED averaged over the 14 deg 2 H-ATLAS SDP 
field and constructed by correlating SPIRE pixel intensities with 100 /im 
IRAS pixel intensity to measure the relative sub-mm colors at each of the 
wavelengths (see text for details) and then scaling the IRAS mean flux over 
the field to SPIRE bands. The best-fit dust temperature and spectral emis- 
sivity parameters are 7^ = 19.0 ± 2.4 K and /} = 1.4 ± 0.4, respectively 
(black solid line). For reference, we also show the expected SED following 
Schlegel et al. (1998) dust map and Finkbeiner et al. (1999) frequency scal- 
ing that involves two temperature components (see Section 4.5; blue solid 
line). 



smaller than the errors quoted above. Figure 2 shows points which 
show zero fluxes in SPIRE maps or are significantly negative; these 
are associated with a combination of the source mask and pixels 
that were associated with parts of the time-streams that are either 
contaminated or contain glitches that were removed. The correla- 
tion analysis above accounts for such pixels when degrading the 
resolution and we find that such points do not bias the color and 
zero-point offset estimates we have quoted above. 

The mean intensity of the SPIRE maps in the SDP field, once 
corrected for So, is around (4.8,2.5, 1.6) (±0.5) MJy sr 1 at 250, 
350, and 500 jim, respectively. These can be compared to the esti- 
mated extragalactic background intensity at each of the three fre- 
quencies of 0.85, 0.69, and 0.39 MJy sr~' with an uncertainty at the 
level of 0.1 MJy sr 1 (Fixsen et al. 1998). The SDP field, on aver- 
age, is a factor of 4 to 5 brighter than the extragalactic background. 
The lowest pixel intensity values of the SDP field in Figure 2 (once 
corrected for So) allow an independent constraint on the extragalac- 
tic background intensity, but such a study is best attempted in fields 
where the overall cirrus intensity is similar to or smaller than the 
expected extragalactic background. Fields that span over a wide 
range of Galactic longitudes and latitudes are also desirable since 
such fields allow an additional constraint on determining a constant 
intensity that is independent of the location. We will attempt such 
studies in future works making use of multiple fields in H-ATLAS. 

Beyond the mean intensity, the extragalactic background aris- 
ing from sources below the confusion noise has been shown to 
fluctuate at the few percent level at 30 arcminute angular scales 
(Amblard et al. 2010). Those faint sources are also responsible for 
roughly 85% of the extragalactic background intensity (Clements et 
al. 2010; Oliver et al. 2010). We are not able to account for the con- 
tamination coming from the unresolved extragalactic background 
light, but the fluctuation intensity of 0.1 MJy sr~' at 30 arcminute 
angular scales do not bias the measurements we report here. The 
background fluctuations act as an extra source of uncertainty in 
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our measurements as they introduce an extra scatter in the inten- 
sity measurements from one region to another and that scatter is 
captured in the overall error budget. 

4.2 Average Dust SED 

Combining the IRAS 100 fim average cirrus intensity over the 
whole SDP area (1.77 MJy sr _1 ) with the above sub-mm color fac- 
tors, and assuming a 15% flux uncertainty coming from the overall 
calibration of SPIRE (Swinyard et al. 2010), we estimate the dust 
temperature T A and /? to be 19.0 ± 2.4 K and 1.4 ± 0.4, respectively 
(Figure 3). The same SED can also be described by the two tem- 
perature model of Finkbeiner et al. (1999), where we fix T d = 9.2 
and 16.2 K and /? = 1.67 and 2.7, with the SED computed using 
the maps derived from Schlegel et al. (1998) dust map following 
an analysis similar to the above (see Section 4.5 for more details). 
The overall fit, however, has a reduced x 1 value of 1 .4 compared 
to 0.8 for the isothermal case. The Polaris flare studied in Miville- 
Deschenes et al. (2010) has a mean intensity of 40 MJy sr~' at 250 
/Jin and 7d = 14.5 ± 1.6 K, /3 = 2.3 ± 0.6 when averaged over the 
flare. The SDP field of H-ATLAS is at least a factor of 8 fainter in 
the mean intensity and has a higher temperature, but a lower value 
for /3 suggesting that the dust temperature and [} vary substantially 
across the sky depending on the intensity of the dust. This compli- 
cates simple approaches to Galactic dust modeling with one or two 
temperatures and /3 values (e.g., Finkbeiner et al. 1999). We discuss 
this further in Section 4.5. 



4.3 SEDs of Bright Cirrus Regions 

Instead of the average (T^,^) values over the whole field, we also 
study the SEDs in five bright cirrus regions that are identified 1 
to 5 in Figure[TJ The previous SED measurement made use of the 
SPIRE maps that were corrected for sub-mm color, G, and So val- 
ues obtained by correlating with the IRAS 100 yum map over the 
whole field and at the pixel scale of IRAS. To study the SED of 
bright regions, we now consider a differential measurement so that 
we can determine the SED independent of the IRAS 100 inten- 
sity in the SDP field. To account for both the extragalactic back- 
ground and the zero-point, we remove from the three SPIRE maps 
the mean intensity from the two low-cirrus regions identified with 
rectangles in Figure 1 . These two regions have intensities that are 
at the low end of SPIRE intensities plotted in Figure 2 and, with 
the previous zero-level included, these intensities are 1.6, 1.2, and 
0.9 (±0.1) MJy sr 1 at 250, 350, and 500 fim, respectively. We do 
the same for the IRAS map and remove the mean intensity of 1 .9 (± 
0.1) MJy sr~' at 100 /im determined for the same two regions. This 
is necessary to avoid introducing an unnecessary difference in the 
relative calibration between SPIRE and IRAS. We account for the 
uncertainty in this mean removal in our overall error budget. 

This procedure allows us to treat the IRAS intensity indepen- 
dent of SPIRE, but the results we show here do not strongly depend 
on this additional step. When we simply used the IRAS S o corrected 
maps, with SPIRE zero-level fixed to IRAS, we still recover SEDs 
that have /? and r d values consistent within uncertainties. The new 
maps, however, lead to differences in the amplitude due to the over- 
all shift in the intensity scale. Another way to think about this is 
that /3 and T d estimates extracted from the isothermal SED depend 
on the intensity ratios and not the absolute intensity. The noise is 
computed by averaging the noise intensity of pixels in the same 
region as defined for the intensity measurements. We add quadrati- 
cally the flux error, the error in the intensity removed from the two 
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Figure 4. Top: Best-fit values of the spectral emissivity parameter, 8, and 
dust temperature, Tj, over the whole SDP field with all maps repixelized at 
6' pixels. The error bars show the l<x uncertainties of these two parameters 
in each of the pixels. The dashed lines show the best-fit models for T L \ vs. 
/3 based on the two analytical descriptions outlined in Section 4.4 with pa- 
rameter values from Dupac et al. (2003) and Desert et al. (2008). The solid 
lines show the same model descriptions but with the parameters obtained in 
this work. We also compare the best-fit line of Veneziani et al. (2010) from 
observations of high density clouds in the BOOMERanG-2003 CMB field. 
Compared to previous measurements that studied high-density dust regions, 
we find a stronger anti-correlation for diffuse emission. The noise-weighted 
mean values of T<$ and 8 are 20. 1 + 1.6K and 1.5 ± 0.3, respectively. These 
are somewhat better than the previous estimates over the whole area as we 
exclude pixels which have relative errors for and /3 that are greater than 
30%. The two points show T& and 8 of Finkbeiner et al. (1999) model-8 
used to scale the dust map of Schlegel et al. (1998). Bottom: Comparison 
between best-fit T<\ and ji distribution with same parameters in the same 
pixels from the Schlegel et al. (1998) dust model (see text for details). 



low-cirrus regions, and an overall calibration error taken to be 15% 
of the intensity (Swinyard et al. 2010). The intensity values of the 
five selected regions are summarized in Table 1. We find the dust 
temperature of these bright cirrus regions are around 18 ± 3K with 
p around 1.5 ± 0.5, and consistent with results found for large-scale 
cirrus observations at high Galactic latitudes with a dust tempera- 
ture of 17.5K (Boulanger et al. 1996). 
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4.4 Relation between fi and T d 

Beyond the bright regions, we also establish the dust temperature 
and modified spectral emissivity parameter over the whole SDP 
area. This allows us to produce maps of Tj and fi over the SDP 
field. To this aim we repixelize all maps to 6' x 6' pixels result- 
ing in a grid of 55 by 55 pixels for the H-ATLAS SDP field; with 
smaller pixel sizes we get noisier estimates of 7^ and fi in regions of 
low cirrus intensity, while with larger pixels we do not sample the 
map adequately. The reprojection is done such that the flux is pre- 
served when going from smaller pixels to 6' pixels. This particular 
choice of pixel size was made so that the total number of grid points 
used for (A, T^,fi) model fits can be completed in a reasonable time 
(a few days in this case) given the computational costs associated 
with generating separate MCMC chains. Here again we consider a 
differential measurement and remove the mean intensity estimated 
in low-cirrus parts (before the maps were repixelized) and in those 
regions we set A = with the assumption of no cirrus. 

With this analysis we construct the two maps that we show 
in Figure 5, where we subselect 6' pixels where the SED fits gave 
relative errors less than 30% for both parameters fi and T^. These 
values are also plotted in Figure 4. They mostly span the high inten- 
sity region of the SDP map. For reference, In Figure 8 (left panels) 
we show the maps of T$ and fi without this selection imposed. At 
first glance fi and Tj appear to be negatively correlated. Part of this 
correlation results from the functional form of the Equation [3] that 
we used to fit the intensity data, especially in the presence of noise 
(Shetty et al. 2009a). 

In order to discriminate any physical effect from the analytical 
correlation, we test two possible models that describe the fi versus 
T d dependence, either with 

P(T A ) = NT^ (4) 
following Desert et al. (2008) or 

PiTi) = — -— , (5) 
C + xl d 

from Dupac et al. (2003). 

Instead of simply using the best-fit and fi values and their 
variances when doing a fit to the two forms of p{T A ) given above, 
we need a numerical method that also takes into account the full 
covariance between the two parameters at each of the pixels. To 
achieve this we fit for the two parameters describing each of the 
relations between fj and T A by making use of the full probability 
distributions captured by the MCMC chains. The procedure we use 
is the same as that of Veneziani et al. (2010). 

A basic summary of the approach is that we fit, for example 
N and a, by drawing random pairs of fi and T d by sampling their 
likelihood functions from the MCMC chains we had first generated 
by fitting the isothermal SED models to individual pixel intensities; 
here again, we restrict the analysis to chains where fi and T4 are de- 
termined with relative errors better than 30% for both parameters. 

By using the full MCMC chains to sample fi and directly 
we keep information related to the full covariance and this takes 
into account the fact that fi and T A are anti-correlated in each of the 
pixels that we use for this analysis. For each of the two forms of 
y8(T d ), we sample the chains by drawing 20,000 random pairs of fi 
and r d ; we established a sampling of 20,000 is adequate by a se- 
ries of simulations using anti-correlated data points in the fi - Ta 
diagram with errors consistent with Figure 4 and assuming random 
correlation coefficients of -0.3 to -0.8. Through this fitting proce- 
dure, we extract the distribution functions of the four parameters N, 
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Figure 5. The dust temperature (top) and spectral emissivity parameter /? 
(bottom) based on our SED fits to pixel intensities that survive the require- 
ment that the two parameters be measured with relative uncertainty (ltr) 
better than 30%. 

a, C, and x and these in return allow us to quote their best-fit values 
and errors. 

We find [N, a] of [N = 1 16 ± 38, a = - 1.4 ± 0. 1] and [C, x] of 
[C = -0.36±0.02, x = (5.1 ±0.1)x 10~ 2 ]. We show these two best- 
fit lines in Figure 4. The two models have^ 2 per degree-of-freedom 
values of 0.8 and 0.5, respectively, suggesting that the model of Du- 
pac et al. (2003) is slightly preferred over the other. With reduced^ 2 
values less than one, it is likely that we are also overestimating our 
overall error budget, especially with the 15% flux calibration un- 
certainty. In Swinyard et al. (2010) the calibration error for SPIRE 
data is stated with an additional safety margin and the likely error is 
between 5% and 10%. For comparison, Desert et al. (2009) found 
[N = 1 1.5 + 3.8, a = -0.66±0.05], while Dupac et al. (2003) found 
C = 0.40 ± 0.02 and x = (7.9 ± 0.5) x 10~ 3 . These two lines, as well 
as the best-fit line of Veneziani et al. (2010), are shown in Figure 4 
for comparison. The Veneziani et al. (2010) measurements involve 
7 high density clouds in the BOOMERanG-2003 CMB field and 
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Figure 6. Example SEDs showing SPIRE and IRAS 6' pixel intensities and 
the best-fit isothermal model with T A = 20.1 and/} = 1.3 (top), T d = 18.2 
and/3 = 2.1 (middle), and T d = 13.1 and/3 = 4.0 (bottom). 



their measurements are consistent with prior works, except for one 
cloud with a low dust temperature of (6.5±2.6)K and /3 of 5.1 ± 1.8. 

While we find f} values as high as 4 to 5, most of the values 
are between 1 and 3. We show three example SEDs in Figure 6 
spanning low, mid, and high values of both T A and /?. Even if we 
constrain the study of P(Ta) relation to the range of 1 < /J < 3, 
we still find non-zero values for the four parameters N, a, C, and 
x consistent with above values, suggesting that there is an intrinsic 
anti-correlation and not driven by the few high ft points. 

Figure 4 demonstrates two interesting scientific results: (i) we 
find an underlying relation between T d and /? that cannot be due 
to an anti-correlation induced by noise when fitting the SED form 



to pixel intensities and (ii) we find a higher value of /3 at low dust 
temperatures compared to the values sugested by previous relations 
in the literature. Our result shows that the anti-correlation also ex- 
ists for low-intensity Galactic cirrus and is not limited to high den- 
sity clumps and molecular clouds that were previously studied. It 
is likely that the T A vs. /3 relation captures different physical and 
chemical properties of the dust grains, including the size distribu- 
tion and the interstellar radiation field that heat the dust. There is 
also a possibility that this anti-correlation results from line of sight 
projection of different dust temperature components (Shetty et al. 
2009b). With the field at a Galactic latitude of 30 degrees, such 
overlap is likely to be small, but perhaps not completely negligible. 
Unfortunately we do not have a way to constrain the line-of-sight 
projection due to the lack of distance information. Also, studies on 
the f} and T A relation are so far limited to handful of fields and more 
work is clearly desirable. We note that laboratory measurements 
have suggested the possibility of such an anticorrelation for cer- 
tain types of dust grains (e.g., Agladze et al. 1996; Mennella et al. 
1998; Boudet et al. 2005). This has been explained as due to quan- 
tum physics effects on the amorphous grains, such as due to two- 
phonon processing and tunneling effects between ground states of 
multi-level systems. Whether the relation we have observed is due 
to averaging different values of temperature along the line-of-sight 
or due to an intrinsic property of dust is something that will remain 
uncertain. 

Related to the observation (ii) outlined above, the measure- 
ments we report here are primarily dominated by the diffuse cirrus 
emission over the whole SDP area. The Dupac et al. (2003) re- 
lation was for a large sample of molecular clouds in the Galaxy 
while the Desert et al. (2008) measurements involve a sample of 
cold clumps detected as point sources in the Archeops CMB ex- 
periment. The dust size distribution is expected to be different in 
denser regions compared to that in diffuse cirrus as the small grains 
are expected to coagulate into large aggregates. The diffuse cirrus 
is likely dominated by small grains and this difference could be 
captured in terms of different values of [} for a given T&. The expec- 
tation is that denser regions would show smaller values of /? (e.g., 
Ossenkopf & Henning 1994), consistent with Figure 4. Once Her- 
schel imaging data have been obtained for more of the H-ATLAS 
areas it will be interesting to study the T d and /? relation for a variety 
of source structures, from dense cores and clumps in our Galaxy to 
extragalactic sources to diffuse Galactic cirrus emission in order to 
establish how the relation changes with the environment. Both the 
Hi-GAL survey with Herschel (Molinari et al. 2010) and Planck 
can make important contributions to this topic in the future. 



4.5 Comparison to a dust model 

We compare our fj and Tj maps with analogous maps obtained 
through model 8 of Finkbeiner et al. (1999) following the dust 
model of Schlegel et al. (1998; SFD hereafter). The SFD dust map 
is produced by combining 100 /mi IRAS and 240 /mi COBE/DIRBE 
data and is an all-sky map of sub-millimeter and microwave emis- 
sion of the diffuse interstellar dust. The model 8 of Finkbeiner et 
al. (1999) involves two dust temperature components at 9.2 and 
16.2 K with /? = 1.67 and 2.70, respectively (see the two points 
in Figure 4). The SFD dust map with a frequency scaling such as 
model 8 is heavily utilized by the CMB experimental community 
both in planning and quantifying the Galactic dust contamination 
in CMB anisotropy measurements. At tens of degree angular scales 
and at frequencies above 90 GHz, Galactic dust is expected to be 
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Figure 7. Left: SPIRE 250 f/m map of H-ATLAS SDP field at the 6' pixel resolution used for studying /J and over the whole map (Section 4.4). Center: A 
map of the same field extracted from the Schlegel et al. (1998) dust map with Finkbeiner et al. (1999) model-8 frequency scaling to obtain 250 /jm intensities 
at the same 6' pixel resolution. Right: A map showing the difference between the observed data and the extrapolation from the Schlegel et al. (1998) dust map 
with Finkbeiner et al. (1999) model-8 frequency scaling. All three maps are in units of Jy/sr. To highlight the difference map, the third map has a shorter color 
stretch. 
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ponent primarily impacts the mm-wave data. We find that the mid 
values are safely produced by the SFD model, but is not a reli- 
able description for regions that have either low or high temper- 
atures. Such regions have either low or high ft values due to the 
anti-correlation between the two parameters. Wide-field imaging 
with Herschel, such as the existing H-ATLAS and the proposed 
Herschel-SnRE Legacy Survey (Cooray et al. 2010) and Planck, 
will provide necessary information to improve the Galactic dust 
map and the associated frequency scaling model as a function of 
the sky position. While in this work we only considered the 14 
deg 2 SDP field, the SGP and NGP portions of H-ATLAS each 
cover about 200 deg 2 , will significantly improve the IRAS-based 
dust model of our Galaxy, especially at Galactic latitudes probed 
by CMB experiments. In a future paper we will return to a further 
analysis on the improvements necessary for the dust model with 
data in those wide fields. 



Figure 9. Angular power spectrum of the H-ATLAS field associated with 
Galactic cirrus. The power-law model fits are considered down to k of 0.03 
arcnrkr 1 to avoid the cirrus measurements with extragalactic background 
fluctuations that are expected to peak at k ~ 0.1 arcnrkr 1 . 



the dominant foreground contamination, especially for polarization 
measurements of the CMB (e.g., Dunkley et al. 2009). 

In order to compare the measurements in the SDP field to pre- 
dictions from the SFD map combined with Finkbeiner et al. (1999) 
frequency scaling for dust emission, we make a new set of maps at 
100, 250, 350, and 500 fim using the SFD dust map and overlap- 
ping with the H-ATLAS SDP field (Figure 7 shows a comparison of 
maps at 250 fim). We analyzed these four simulated maps by apply- 
ing the same procedure as we used for extracting spectral emissivity 
parameter and dust temperature with SPIRE and IRAS real maps. 
We also include the SPIRE noise making use of the SPIRE noise 
maps of the field generated from the data. We obtain two maps from 
the SFD dust map, one connected to the spectral emissivity (3 and 
one related to the dust temperature T^. 

We show the ratio map between data and SFD model results 
for both spectral emissivity and dust temperature in Figure 8 (right 
panels), while the left panels show our measurements directly on 
SPIRE data. The two component dust model of Finkbeiner et al. 
(1999) only captures a limited range of dust temperature and p. The 
model involves a low-dust temperature component at 9.2 K with 
P = 1.67; SPIRE maps spanning out to 500 jim are not strongly 
sensitive to such a cold dust component though such a cold com- 



5 CIRRUS POWER SPECTRUM 

Given that SPIRE is capable of mapping the diffuse emission at 
large angular scales, we also study the angular power spectrum re- 
lated to cirrus emission for each of the wavelength bands. For this 
measurement we keep the maps at the original pixel scale (Sec- 
tion |2) and compute the power spectrum of the intensity in maps 
masked for detected sources. With 1% to 5% of the pixels masked, 
we found the mode coupling introduced by the source mask to be 
negligible. 

We show our measurements in Figure 9. At k > 0.1 arcmkr 1 
the effects of the beam transfer function become important but we 
do not make a correction for the beam here as we are mostly in- 
terested in the power spectrum at degree angular scales, where the 
transfer function related to the SPIRE beam is effectively one (Mar- 
tin et al. 2010; Miville-Deschenes et al. 2010; Amblard et al. 2010). 
In addition to the beam, there is also the map-making transfer func- 
tion associated with any filtering employed during the map-making 
process (Pascale et al. 2010). At the angular scales of interest this 
transfer function is also consistent with one. 

We assume a power-spectrum of the form P(k) = Po(k/koY to 
describe the measurements and take k = 0.01 arcmin 1 to be con- 
sistent with previous studies. Using data out to k < 0.03 arcmkr', 
to avoid contamination with fluctuations associated with the ex- 
tragalactic background, we find y = -2.4 ±0.1 at 100 /jm with 
IRAS and -2.6 ± 0.2 at each of the SPIRE bands at 250, 350, 
and 500 /jm. The measured values of the fluctuation amplitudes 
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Figure 8. Top: Left - Map of the dust temperature: best fit values are shown. Units are degrees Kelvin. Pixels with < 5 K are dominated by noise and 
should not be considered as an accurate measurement. Figure 5 shows the subselection where Ti and /? are determined with a relative accuracy better than 
30%. Right - Ratio map between data temperature map and a map of the dust temperature by reanalyzing the Schlegel et al. (1998; SFD) dust map for the 
same region scaled by model-8 of the Finkbeiner et al. (1999) two temperature description and refitted with our isothermal SED following the same procedure 
as data. Bottom : Left - Map of spectral emissivity parameters /?, best fit values are shown. Right - Map of the ratio between data spectral emissivity map and 
the pectral emissivity map derived from the SFD model dust map similar to dust temperature. 



are P = (8.3 ± 0.7) x 10 6 , (2.2 ± 0.2)xl0 7 , (4.8 ± 0.5) x 10 6 , 
( 1 .3 ± 0. 1) x 10 6 at 100, 250, 350, and 500 fim respectively (in units 
of Jy 2 /sr). The normalization at 100 arcmkr 1 angular scale we find 
for IRAS 100 yum map is fully consistent with the relation between 
P(k = 0.01) and the mean 100 jim intensity in Miville-Deschenes et 
al. (2007; see their Figure 4), given the IRAS mean intensity of 1.77 
MJy sr~' in the SDP field. The power-law slope y we find is some- 
what lower than measurements for the slope in the literature (e.g., 
Miville-Deschenes et al. 2010 with a slope if -2.7 ± 0.1, Martin et 
al. 2010 with slopes of -2.74 ± 0.03 and -2.81 ± 0.03 in two sepa- 
rate fields), but consistent with the analysis in Roy et al. 2010 (with 
a slope of -2.6 ± 0.07 at 250 yum). All these measurements are con- 
sistent with each other given the overall uncertainties. It could also 
be that our slope is lowered by tens of percent level due to fluctua- 



tions associated with the extragalactic background. If we constrain 
to k < 0.01 arcmin 1 , keeping only two data points, we do find a 
higher slope closer to -2.8 to -2.9 but with larger uncertainty (±0.4). 

The -\/k 2 P(k)/2n captures the rms fluctuations over the whole 
SDP area arising from Galactic cirrus and the extragalactic back- 
ground at a given value of the wavenumber. At large angular scales, 
the fluctuations generated by the extragalactic sources (mainly the 
sources contributing the background confusion noise) are subdom- 
inant with values of order 10 4 to 10 s Jy 2 /sr in P (k) when k < 0.01 
arcmin -1 (say at 250 (im; Amblard et al. 2010). For comparison, in 
Figure 9 we find Galactic cirrus fluctuations at the level of 10 7 to 
10' Jy 2 /sr. Thus, at k = 0.01 arcmin 1 we can safely assume that 
all of the fluctuations we have measured arise from Galactic cirrus. 
Taking the VCo) values we can make an independent estimate of 
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the dust temperature and p similar to the analysis of Martin et al. 
(2010). We find T A = 20.1 ± 0.9 K and = 1.3 ± 0.2, consistent 
with the same two quantities we obtained in Section 4.2 by cross- 
correlating SPIRE maps with IRAS. 



6 SUMMARY AND CONCLUSIONS 

In this paper, we have studied the Galactic dust SED and the an- 
gular power spectrum of dust fluctuations in the 14 deg 2 Science 
Demonstration Phase (SDP) field of H-ATLAS. By correlating the 
SPIRE 250, 350 and 500 fim and IRAS 100 fjm maps to extract the 
sub-mm color terms of SPIRE maps relative to the IRAS 100 yum 
map of the SDP field, we find the average dust temperature over 
the whole field to be 19.0 ± 2.4 K with the spectral emissivity pa- 
rameter taking a value of 1.4 ± 0.4. We find all of the bright cirrus 
regions to have dust temperatures T d over a narrow range of 17.4 
to 18.3 K (±2.5 K), with a spectral emissivity parameter ft ranging 
from 1.4 to 1.9 (±0.5). Similar to previous studies, we find an anti- 
correlation between T A and P; when described by a power-law with 
P = AT a d , we find A = 116 ± 38 and a = -1.4 ± 0.1 while a 
relation of the form j3 = (C + xT A y l is also consistent with data 
with C = -0.36 ± 0.02 and x = (5.1 ± 0.1) x 10~ 2 . The observed 
inverse relation between T A and p is stronger than the previous sug- 
gestions in the literature and we have suggested the possibility that 
this stronger anti-correlation may be due to the fact that we study 
primarily diffuse cirrus while previous studies involved high den- 
sity environments such as molecular clouds and cold clumps. We 
also make an independent estimate of the dust temperature and the 
spectral emissivity parameter, when averaged over the whole field, 
through the frequency scaling of the rms amplitude of dust fluctu- 
ation power spectrum. At 100 arcminute angular scales, we obtain 
T A = 20.1 ± 0.9 K and /? = 1.3 ± 0.2, consistent with previous esti- 
mates. The cirrus fluctuations power spectrum is consistent with a 
power-law at 100, 250, 350 and 500 fim with a power-law spectral 
index of -2.6 ± 0.2 from 1 to 200 arcminute angular scales. 

After we completed this paper, we became aware of a simi- 
lar study involving the p and T$ in the first two fields covered by 
the Hi-GAL survey (Paradis et al. 2010). These authors also find 
an anti-correlation between the two parameters, but the mean re- 
lation is distinctively different between the two fields at Galactic 
longitudes of 30 and 50 degrees, with both on the Galactic plane. 
When compared to the H-ATLAS SDP field, these two fields have 
mean intensities that are a factor of 200 larger at the level of 1000 
MJy sr~'. Interestingly the P(T A ) relation they find for the field at 
I = 30 degrees is consistent with the relation we report here when 
extrapolating their relation that was determined over the range of 
1.5 < P < 2.5 and 18 < T d /K < 25 to lower 14 K dust tempera- 
tures we find in some of the pixels in our field with /? ~ 4. It could 
very well be that the dust properties are far more complex than the 
simple isothermal models we have considered and a variety of ef- 
fects may be contributing to the observed anti-correlation. Further 
studies making use of wide area maps are clearly warranted. 
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